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Osl ■ ABSTRACT 

, Context. For the past three decades, and until recently, there has been a serious discrepancy between the observed and theoretical 

r-y ■ values of the apsidal motion rate co of the eccentric eclipsing binary DI Her, which has even been interpreted occasionally as a 

' possible failure of General Relativity (GR). A number of plausible explanations have been put forward. Recent observations of the 

1/^ ' Rossiter-McLaughlin effect have shown convincingly that the reason for the anomaly is that the rotational axes of the stars and the 

orbital axis are misaligned, which changes the predicted rate of precession significantly. 

Aims. Although as a result of those measurements the disagreement is now drastically smaller, it remains formally at the level of 50%, 
possibly due to errors in the measured apsidal motion rate, outdated stellar models, or inaccuracies in the stellar parameters. The aim 
' of this paper is to address each of these issues in order to improve the agreement further. 

, Methods. New times of minimum have been collected, and used for a redetermination of the apsidal motion rate. Based on the 

latest determinations of the absolute dimensions of the binary, we have computed new stellar evolution models with updated physical 
inputs, and derived improved apsidal motion constants for the components. We have performed Monte Carlo simulations to infer the 
S-^' theoretical distribution of io, including the contributions from GR as well as tidal and rotational distortions. All observational errors 

Q , have been accounted for. 

5_j ■ Results. Our simulations yield a retrograde apsidal motion rate due to the rotationally-induced oblateness of -0.00056 deg cycle"' 

' (mode of the distribution), a GR contribution of +0.00068 deg cycle"', and a tidal contribution of +0.00034 deg cycle"', leading to a 

^ ' total predicted rate of +0.00046 deg cycle"'. This is in excellent agreement with the newly measured value of +0.00042 deg cycle"'. 

The formal difference is now reduced to 10%, a small fraction of the observational uncertainties. 

Key words, stars: binaries: close; stars:evolution; stars:interiors; stars: fundamental parameters; stars:rotation; General Relativity 

> 

^ ■ 

1. Introduction 2. The presence of a circumstellar cloud between the compo- 

nents, which would affect the gravitational field of the stars; 
(N ; For some three decades the eccentric B -type eclipsing binary 3 ^.^^^^^^ ^ ^.^^^^^ ^^^^ ^ ^^.^^ 

^ . DI Her (P = 10.55 days) has been known to have an apsi- ^^^j^ contribute to the motion of the line of apsides; 

Q , dal motion rate that is too slow compared to theoretical predic- 4 alternative theory of gravitation. 

tions (for a brief review, see Claret 1998). Until very recently 

^ , this presented a serious challenge to our understanding of stel- Some of these ideas were examined quantitatively and ruled out 

— ■ lar physics, and even the vaUdity of General Relativity has been ^y Claret (1998). The third body hypothesis has been a serious 

called into question given that the GR contribution is domi- contender for many years (see, e.g., Khodykin 2007), although 

nant for DI Her. The discrepancy in the total rate of precession spectroscopic or photometric evidence of it has ever been 

jj amounts to 400% or more. This anomaly persisted even when found. Difficulties in the measurement of such a small apsidal 

cq ^ considering hypothetical stellar models with infinite mass con- motion ^s that of DI Her (wobs = 1-04° per 100 yr, or 0.00030 

centration (apsidal motion constant ^2 = 0), which tends to re- jeg cycle"' ; Guinan et al. 1994) were also suggested as a possi- 

duce the predictions, and this was often interpreted as a failure of ^jg explanation from time to time 

General Relativity. The problem was all the more puzzling when ^n atti-active alternative is a configuration for the binary in 
considering that studies of other binary systems suitable for test- ^^ich the spin axes are tilted relative to the orbital axis (Shakura 
ing GR and with well determined absolute dimensions indicated 1985^ Q^jn^n & Maloney 1985, Company et al. 1988, Claret 
good agreement with theory (see, e.g., Claret 1997 and Wolf et 1993) ^his would introduce a correction to the rotational term 
al. 2009). This suggested to many that the problem may he with j^e apsidal motion that can be negative (regression of the line 
DI Her itself. A number of hypotheses were put forward over the apsides), resulting in a smaller total co. The young age of the 
years in order to explam the disagreement, including system (a few Myr) certainly makes this plausible, as tidal forces 
1. A rapid ongoing circularization of the orbit, affecting the ^o^l^ have had time to operate and bring the spin axes of 
measured times of minimum- ^^^^ alignment with the axis of the orbit (see Sect. |3Jl. 
' A remarkable confirmation of this prediction was obtained re- 
Send offprint requests to: A. Claret, e-mail:claret@iaa.es cently by Albrecht et al. (2009), who provided compelling obser- 
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vational evidence for a peculiar orientation of the stars by taking 
advantage of the Rossiter-McLaughlin effect. This is a distortion 
of the spectral lines that occurs during the eclipses, as a result 
of the partial obscuration of the rotating stellar disks. Using this 
effect, they were able to measure the angles between the sky pro- 
jections of the spin axes and the orbital axis, and show that there 
is indeed a strong misalignment for both components of the bi- 
nary. It is the first documented case in which this has been shown 
convincingly for stars in a double-lined eclipsing binaryQ 

With this new information Albrecht et al. (2009) managed to 
reduce the discrepancy with the predicted apsidal motion rate of 
DI Her very significantly, so that theory and observation are now 
in much better agreement. Some of the angles involved in com- 
puting the effect of the tilted axes remain unknown, however, so 
the theoretical prediction necessarily assumes a statistical char- 
acter in this case, which Albrecht et al. (2009) tackled by means 
of Monte Carlo simulations. The nominal difference between the 
measured apsidal motion rate and the mode of their theoretical 
probability distribution for io is still about 50% at present, al- 
though the distribution is rather wide so theory and observation 
are formally consistent with each other. Some of this difference 
could be due to the measurement of ioobs itself, to inaccuracies 
in the adopted absolute dimensions of the stars, or to the stellar 
evolution models used to infer the apsidal motion constants k2, 
which are needed as inputs to the theoretical calculations. 

In this paper we re-examine all the ingredients in order to 
provide the most accurate predictions and observational con- 
straints. We incorporate new measurements of the times of 
eclipse to derive an improved value of titobs for DI Her that is sig- 
nificantly larger than the value adopted by Albrecht et al. (2009). 
We use the most recent determinations of the mass, radius, and 
other properties of the stars, and we revisit the temperature de- 
terminations. With these we provide new fits to stellar evolution 
models and derive new apsidal motion constants. We then repeat 
the comparison of the measured apsidal motion rate with theory, 
also examining the sensitivity of the predictions to several of the 
key observables. 



2. Observational data 

2. 1 . The absolute dimensions 

The absolute masses and radii of the components of DI Her used 
here are adopted from the recent compilation by Torres et al. 
(2009a), in which these and other properties were re-examined 
or redetermined with a uniform methodology. The uncertain- 
ties are adopted from the same source. We summarize all ob- 
servational results in Table [1] and describe some of them later. 
The effective temperatures reported by Torres et al. (2009a) are 
based on the individual B-V colors as given by Popper (1982), 
and on his color/temperature relation (Popper 1980). As a check, 
we used the Stromgren photometry obtained by Hilditch & Hill 
(1975) and the calibration of Crawford (1978), and derived a 
mean de-reddened color for the binary that is essentially identi- 
cal to Popper's. The color/temperature calibration by Paunzen 
et al. (2005), which combines several more recent tempera- 
ture scales, then yields individual temperatures of 17300 K and 
15400 K for the primary and secondary, with uncertainties of 
800 K. These are the temperatures we adopt in this work, which 
are 300 K hotter than those listed by Torres et al. (2009a). 



' Spin-orbit misalignment has been reported previously for several 
host stars of extrasolar transiting planets. 



Table 1. Observational parameters for DI Her. 



Parameter 


Primary 


Secondary 


Mass m (Mq) 


5.17 + 0.11 


4.524 ± 0.066 


Radius R (Rq) 


2.681 +0.046 


2.478 ± 0.046 


Effective temperature (K) 


17300 ± 800 


15400 ± 800 


V sinyS (km s"') 


108 + 4 


116 + 4 


A (deg) 


+72 ±4 


-84 ±8 


logfc 


-2.10 ±0.05 


-2.13 + 0.05 


Orbital period P (days)* 


10.5501696 ± 0.0000007 


Reference epoch Tq (HJD) 


2,447,372.9567 


+ 0.0005 


clio at reference epoch (deg) 


330.0 


+ 0.1 


Orbital eccentricity e 


0.4895 ± 0.0008 


Orbital inclination i (deg) 


89.30 ± 0.07 


diobs (deg cycle"') 


0.00042 + 0.00012 



* P corresponds here to the sidereal period P,. 



2.2. New apsidal motion rate 

The apsidal motion rate of DI Her was determined by means 
of an 0-C diagram analysis. The method we followed is that 
described by Gimenez & Garcia-Pelayo (1983), which is an it- 
erative weighted least-squares procedure including terms in the 
eccentricity up to the fifth order There are five adjustable quanti- 
ties to be determined: the rate of periastron advance d>, the side- 
real period f j, the eccentricity e, and the longitude of periastron 
ajQ at the reference epoch Tq. The periastron position ai at any 
epoch E is then given by the linear equation 

tJ — tJQ + d> E , (1) 

and the relation between the sidereal and the anomaUstic periods, 
Ps and P„, is given by 

P,^Pa(l-cj/360°). (2) 

The value for the orbital inclination was adopted as / = 89°30 
from the light-curve analysis of Guinan & Maloney (1985). 

We used all precise photoelectric and CCD timings recently 
recalculated by Kozyreva & Bagaev (2009, their Tables 1 and 2), 
as well as new times of eclipse obtained by us at two observato- 
ries, with the following instrumentation: 

• Observatory and Planetarium of Johann Palisa, Ostrava, 
Czech Republic: 0.20-m or 0.30-m telescopes with an 
SBIG ST-8XME CCD camera and Johnson V, R, and / fil- 
ters; 

• Valasske Mezinci Observatory, Czech Republic: 0.30-m 
Celestron Ultima telescope with an SBIG ST-7 CCD cam- 
era and Johnson R filter 

The reduction of all images was performed with the synthetic 
aperture photometry software developed by Motl (2007). The 
frames were dark-subtracted and flat-fielded, and the heliocen- 
tric correction was applied. An example of these observations 
is presented in Figure [T] The new times of primary and sec- 
ondary eclipse and their errors were determined by the classi- 
cal Kwee & van Woerden (1956) algorithm, and are listed in 
Table |2] Earlier visual and photographic times were not used in 
our analysis because of the large scatter of those measurements. 
A total of 69 times of minimum light were employed, including 
34 corresponding to the secondary eclipse. 

The solution was performed by minimizing the standard 
sum, and the errors of the observations were adopted as pub- 
lished. The resulting apsidal motion parameters and their uncer- 
tainties from the least squares fitting are presented at the bottom 
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Fig. 1. Differential light curve of DI Her during the primary min- 
imum of 2005 May 20, as measured by Ms. Hana Kucakova 
at the Johann Palisa Observatory, Technical University Ostrava, 
Czech Republic. The observations were made through a Johnson 
R filter with a 20-cm reflecting telescope and SBIG ST-8XME 
CCD camera, with exposure times of 25 seconds. The compar- 
ison star was GSC 02109-00247 {V = 10.04). The light curve 
contains more than 500 points obtained during nearly 6 hours of 
monitoring. 

Table 2. New times of minimum light of DI Her. 




1000 



Fig. 2. The 0-C diagram of DI Her (eclipse timing residuals), 
together with our best-fit apsidal motion model. The solid line 
and circles correspond to the primary, and the dashed line and 
triangles to the secondary. 
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-0.003 




HJD-2,400,000 


cr (days) 


Epoch 


Filter 


Observatory 


52456.46107 


0.0005 


482.0 


R 


Ostrava 


53511.47699 


0.0004 


582.0 


R 


Ostrava 


53933.48379 


0.0005 


622.0 


I 


Valasske Mezincf 


54239.43802 


0.0008 


651.0 


R 


Ostrava 


54627.36687* 


0.0005 


687.5 


VRI 


Ostrava 



* Mean value of VRI measurements. 



of Table [T] and the corresponding 0-C diagram is shown in 
Figure|2] The value for the best fit was - 500. This is con- 
siderably larger than the number of degrees of freedom (v = 64), 
which is most likely due to the formal errors of the eclipse tim- 
ings being underestimated, as is often found to be the case. The 
uncertainty of (1) reported in Table [1] already accounts for this ef- 
fect. The extra scatter compared to the uncertainties is seen more 
clearly in Figure |3] in which the influence of the apsidal motion 
terms has been subtracted from the measurements. There is no 
evidence in this diagram of any other significant variations, such 
as might be induced by the presence of a third body in the sys- 
tem. 



3. Stellar evolution models 

As a prerequisite for investigating the theoretical apsidal motion 
of DI Her, the observed stellar properties such as the mass, ra- 
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Fig. 3. The 0-C residuals from the photoelectric and CCD tim- 
ings of DI Her, after subtraction of the apsidal motion terms. 
Symbols are as in Figure|2] 



dius, and temperature of each component must be well matched 
by stellar evolution models, at a single age for the binary, since 
these models are then used to infer the apsidal motion constants 
k2 that factor into the rotational and tidal contributions to (1>. The 
basic properties of the stellar evolution code used here have been 
described by Claret (2004). The models include overshooting be- 
yond the formal boundary of the convective core, as determined 
by the Schwarzschild criterion. The distance over which this is 
considered to occur is specified as dover = ffovH^, where H^ is the 
pressure scale height taken at the edge of the convective core, as 
given by Schwarzschild's criterion, and cKov is a free parameter 
We adopt here q-qv = 0.20, following the results of Claret (2007). 
The evolution code incorporates radiative opacities from the ta- 
bles provided by Iglesias & Rogers (1996), and the calculations 
of Alexander & Ferguson (1994) for lower temperatures. 

The chemical composition of double-lined eclipsing binaries 
is often difficult to determine observationally, and that is the case 
also for DI Her This introduces some uncertainty in the com- 
parison with stellar evolution models, since the metal abundance 
Z then becomes a free parameter. For young systems such as 
DI Her, Z is quite sensitive to the effective temperatures of the 
stars. We computed evolutionary tracks for the measured masses 
and explored a range of compositions, finding the best fit for {X, 
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Fig. 4. Stellar evolution models compared with the measure- 
ments for DI Her. The evolutionary tracks are computed for the 
measured masses (see Table [TJ. The best fit has {X, Z) - (0.71, 
0.02), implying a chemical composition near solar. 
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Fig. 5. Radius evolution of DI Her, based on the models shown 
in Figure |4] The shaded areas represent the measured radii and 
their uncertainties. The best-fit age of about 4.5 Myr is indicated 
by the vertical dotted line. 



Z) - (0.71, 0.02), which suggests an abundance very near that of 
the Sun. The observations along with our best fit are displayed in 
Figure |4] The agreement with the temperatures is quite satisfac- 
tory, given the error bars and other uncertainties described be- 
low. Precise age determinations for young systems are not easy 
to make and require very accurate measurements of the absolute 
radii. This is illustrated clearly in Figure |5] As small as the ra- 
dius errors already are in DI Her (< 2% relative errors, which 
are among the best available for eclipsing binaries; see Torres et 
al. 2009a), the uncertainties allow for a wide range of ages be- 
tween ~2 Myr and ~7 Myr. Nominally the best fit to the radii is 
achieved at about 4.5 Myr with these models. We note, in pass- 
ing, that other models such as those from the Yonsei-Yale series 
(Yi et al. 2001, Demarque et al. 2004) do not yield as good a 
fit, at least with the grid of standard compositions (X, Z) that are 
publicly available. 

The theoretical apsidal motion constants inferred from our 
models are log ;fc2i = -2.10 + 0.05 and log ^22 = -2.13 + 0.05 for 



the primary and secondary, respectively. These values indicate 
that the modern models favor stars that are slightly more cen- 
trally concentrated in mass than those adopted by Claret (1998) 
and Albrecht et al. (2009). This goes in the direction of making 
the predicted apsidal motion rate smaller Tests were made in or- 
der to evaluate the influence on the inferred ^2 of uncertainties 
in the age determination and in the chemical composition, and 
these uncertainties are included in the errors reported above. 

Regarding the measured effective temperatures, it is worth 
mentioning that because the spin axes of the stars are tilted rel- 
ative to the orbital axis, the temperature distribution over their 
disks due to gravity darkening could be highly non-uniform and 
asymmetric as viewed from the Earth, and it is even possible 
one or both components are viewed nearly pole-on. There may 
therefore be a bias in the measured effective temperatures that 
could amount to several hundred degrees for DI Her, but because 
the exact orientation is not known (see next section), the effect 
is difficult to estimate more precisely. Rotation itself changes 
the internal structure of stars and affects their evolution, gener- 
ally leading to cooler temperatures. To explore this effect more 
quantitatively for DI Her we have computed models accounting 
for rotation (rigid-body approximation), and we find differences 
of the order of 200 K, depending on the inclination relative to 
the line of sight. Once again, since the precise orientation of the 
spin axes is unknown, we are unable to be more specific regard- 
ing this effect. 

With the models computed in this section we have calculated 
the age at which the stellar spin axes in DI Her are expected 
to become aligned with the axis of the orbit due to the action 
of tidal forces. The procedure follows closely that described by 
Torres et al. (2009b, and references therein) for the case of the 
eclipsing system LV Her. The result using the theory by Zahn 
(1977, 1989) is ~10^ yr With the more general equations by 
Hut (1981), valid for systems of high eccentricity such as this, 
we obtain 3x10^ yr, which is an order of magnitude larger than 
the evolutionary age. These estimates suggest, as is observed, 
that tidal forces have not had sufficient time in this system to 
make the rotation axes parallel to the axis of the orbit. 

4. The apsidal motion of DI Her 

The total rate of apsidal motion from theory, to be compared with 
the measured value, is given by sum of the individual contribu- 
tions of each component due to tidal and rotational distortions, 
and the general relativistic term. When accounting for the possi- 
bility of misalignment (Shakura 1985), we have 



'i'total = i^GR + Wtidal.l + Wtidal,2 + ^^,-ot,l 01 + Wiot,2 I 

where 

I5m2-j 



•J^tidal,; 



kijgie) 



^) (l + ^^ 

Q.K \ nij 



f(e)k2j 



(3) 



(4) 



(5) 



sin^ / 



cos aj (cos Oj - cos /3j cos /) + - sin^ / (1 - 5 cos^ aj) 



The subindices 'tidal,/ and 'rot,/' in eq.O) denote the tidal 
and rotational contributions of component j, and 'GR' represents 
the relativistic contribution to the periastron advance. The angle 
/ is the measured inclination of the orbital plane relative to the 
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line of sight (see Table[TJ, aj are the angles between the rotation 
axes and the normal to the orbital plane, and /3j are the angles 
between the rotation axes and the line of sight to the observer. 
The symbols Q.j and Q.k represent the rotational angular velocity 
for component j, and the Keplerian orbital angular velocity; A is 
the semimajor axis of the orbit, and Rj and nij are the radius 
and mass of component /'. The symbols k2j represent the apsidal 
motion constant of component j, which we determined in the 
previous section. The auxiliary functions /(e) and g{e) depend 
only on the eccentricity e, and can be written as /(e) = (1 -e^y^ 
and g(e) = (8 + 126^ + e'^)f(ef -^/S. The factors (pj capture the 
dependence on the orientation of the stars. When the spin axes 
are perfectly aligned with the orbit, aj - 0, /3j = /, and (pj = I, 
so that eq.© reduces to the classical expression presented by 
Sterne (1939a, 1939b). 

The relativistic contribution to the apsidal motion, which is 
independent of stellar structure, was given by Levi-Civita (1937) 
as 



WGR = 2.29 X 10" 



,(mi + m2) 



A(l 



deg cycle 



(7) 



where m, and A are in solar units. 

Introducing numerical values in eq.(|4]i, we obtain for the 
tidal conti-ibution Wtidai.i + Wtidai,2 - 0.000339^3° '^^S cycle"', 
where the error bars are given in units of the last significant 
digit and account for all observational uncertainties. This is 
about 10% smaller than obtained by Albrecht et al. (2009), who 
adopted slightly different absolute dimensions for the stars. The 
relativistic contribution is wgr = 0.000677;^^ deg cycle"'. 

The calculation of the rotational contribution is more in- 
volved due to the dependence on the orientation of the sys- 
tem. The angles aj and fij in eq.© are not directly measured. 
However, they are related through the expression 



cos aj - cos f3j cos / + sin /3j sin / cos Aj , 



(8) 



in which Aj are accessible to observation through the Rossiter- 
McLaughlin effect, and represent the angles between the projec- 
tions of the spin axes and the orbital axis on the plane of the 
sky, for each star Albrecht et al. (2009) determined their val- 
ues to be Al = +72° + 4° and A2 = -84° ± 8° for the primary 
and secondary of DI Her, respectively. The angular velocities 
ilj depend on the projected linear rotational velocities, which 
were measured spectroscopically by Albrecht et al. (2009) as 
vi sin/Si = 108 ± 4 km s"' and V2 sinj62 = 116 + 4 km s"', and 
on the unknown inclination angles /3j needed to de-project the 
measured velocities in order to infer the equatorial values. We 
may therefore write 



D.J _ P [vjsm/3j] 



27:Rj sin /3j 



(9) 



in which we represent vj sin /3j in square brackets to indicate that 
it is a measured quantity. 

We follow Albrecht et al. (2009) and perform Monte Carlo 
simulations in order to determine 0i and (f>2, and with these the 
total rate of apsidal motion in eq([3]). The result is a probability 
distribution rather than a single value. For this we assume that 
the angles /3 j are distributed completely at randorrQ, or equiva- 
lently, that cos /3j has a uniform distribution. Knowledge of /3j is 
sufficient to determine cos aj through eq.®, and Q^/Q^: through 
eq.®. As did Albrecht et al. (2009), we exclude angles /3j that 
lead to equatorial rotational velocities exceeding the breakup ve- 
locity, for which we adopt 600 km s"' . 
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Fig. 6. Monte Carlo calculations for the apsidal motion of 
DI Her, based on 2 x 10^ simulations. From top to bottom, the 
first three panels show the separate contributions from GR, tidal 
distortions (2nd and 3rd terms of eq.|[3l), and rotational distor- 
tions (4th and 5th terms of eq.|[3l), based on the measurements 
by Albrecht et al. (2009). The distributions account for all ob- 
servational errors. The bottom panel shows the predicted total 
apsidal motion rate on the left (sum of all the above contribu- 
tions), and on the right the result obtained if the spin axes are 
assumed to be aligned with the orbital axes, rather than strongly 
misaligned, as observed. The measured rate of apsidal motion 
of Wobs = 0.00042 ± 0.00012 deg cycle"' is indicated with the 
vertical dashed line and shaded area, and agrees remarkably well 
with the calculation for misaligned axes. The mode of the distri- 
bution for the case of aligned axes, on the other hand, is more 
than 5 times larger than the measured value. 



In the first three panels of Figure |6] we show separately the 
GR, tidal, and rotational contributions to the apsidal motion 
rate. The distributions for the first two were also generated with 
Monte Carlo simulations in which we have assumed all obser- 
vational quantities have normally distributed errors. The distri- 
bution for the rotational term in the third panel is much wider 
and asymmetrical, but also sharply peaked, and has a negative 
modal value of ciiiot = -0.00056 deg cycle"', as proposed by 
Shakura (1985) and others, indicating a retrograde motion of the 
line of apsides due to the misalignment of the spin and orbital 
axes. This result is numerically some 10% larger than that ob- 
tained by Albrecht et al. (2009), on account of the slightly dif- 
ferent absolute dimensions and theoretical k2j values adopted. 
The bottom panel of Figure |6] displays the predicted distribu- 
tion for the total apsidal motion rate (left), which has a mode of 
Wtotai - 0.000462^^^3 deg cycle"' (asymmetric Icr errors enclos- 
ing 68.3% of the results). The agreement with the new measured 
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value of Wobs = 0.00042 + 0.00012 deg cycle"' from Sect. O 
is excellent, although the theoretical distribution has a long tail 
toward the left that could certainly accomodate smaller values 
of Wobs- Still, it is interesting to note that the measurement now 
agrees much more closely with the most probable value from 
theory, the formal difference being only 10%. The improvement 
compared to the results by Albrecht et al. (2009) is due mostly 
to the refinement of the measured apsidal motion rate, and to 
some extent also to the updated absolute dimensions and stel- 
lar models leading to more realistic apsidal motion constants. In 
particular, our models imply stars that are somewhat more cen- 
trally concentrated in mass. 

If the spin axes of DI Her are assumed to be perfectly aligned 
with the orbit, instead of strongly misaligned as observed, the 
resulting Monte Carlo distribution for the total apsidal motion 
rate is that shown on the right-hand side of the bottom panel of 
Figure|6l The mode of this distribution is at Wtotai - 0.00220^[3 
deg cycle"', which is 5.2 times larger than the measured value. 
The distribution is completely inconsistent with the observation. 
This result for the case of aligned axes is nearly a factor of two 
larger than that presented by Albrecht et al. (2009). The reason 
for this significant difference is that theirs is based on calcula- 
tions by Guinan & Maloney (1985) in which the projected ro- 
tational velocities of both stars were taken to be 45 km s"', as 
measured spectroscopically in the 1970s, whereas we used the 
values reported by Albrecht et al. (2009) themselves, obtained 
in 2008, which are considerably larger (Table [1]). 

The sensitivity of our Monte Carlo calculations for the mis- 
aligned case to some of the quantities measured by Albrecht et 
al. (2009), such as Aj and vj sin /3j, is illustrated in Figure |7] We 
show the changes in the resulting distributions when these val- 
ues are altered one at a time for the primary star, with all other 
observational quantities held at their nominal values. 

The DI Her observations of Albrecht et al. (2009) exploiting 
the Rossiter-McLaughlin effect represent a breakthrough in our 
understanding of this system, and the accuracy of their measure- 
ments has played a central role in achieving the now excellent 
agreement with the measured apsidal motion rate of the binary, 
to within 10% (a small fraction of the measurement eiTor). The 
main factor leading to this improvement has been the new mea- 
surement of the apsidal motion rate presented here, followed by 
other refinements including the consideration of stellar models 
with updated input physics, and more accurate absolute dimen- 
sions for the binary components. We now know that the spin 
axes of the DI Her components are strongly misaligned, and the 
mystery of its slow apsidal motion rate is solved. However, ques- 
tions remain. During the process of binary star formation the 
stellar spins and the orbit presumably derive their angular mo- 
mentum from the same source, which is the parent molecular 
cloud. One would then naturally expect that the axes would be 
closely aligned, yet in this system they are not. Thus, interest 
in this object is sure to continue as theorists try to answer these 
questions about DI Her, and observers look for other systems 
like it. 
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Fig. 7. Sensitivity of the Monte Carlo calculations for the pre- 
dicted apsidal motion rate of DI Her to changes in some of the 
input quantities measured with the Rossiter-McLaughlin effect. 
The hatched histogram is the result using the values exactly as 
reported by Albrecht et al. (2009), and the vertical dashed line 
and shaded region marks the new measured apsidal motion rate 
reported in this work. The shaded histogram on the left is for 
vi sinySi - 160 km s"' instead of the nominal value of 108 
km s"', with all other quantities left at their measured values. 
The shaded histogram on the right is for Ai - 30° instead of 
the observed value of 72°. Similar effects are seen when varying 
those parameters for the secondary. 



